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Abstract 

It is shown that the contribution to the induction which at an internal point of a spin 
density distribution is mathematically described as a local is virtually caused by the 
summing-up of the fields created by all elements of this distribution. Therefore, the 
proportionality coefficient between this contact (Fermi) field and magnetic moment 
density at the point of observation is equal to 8ir/3 only for spherically symmetrical 
s-shells. If the symmetry of spin density distribution lowers, the value of this coef- 
ficient becomes dependent on the spin direction. As a sequence, in low-symmetry 
crystals and molecules additional anisotropic contributions to the hyperfine field 
emerge. PACS: 76.60.Jx 
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Introduction 



Spin magnetic moments of electrons from proper shells of an isolated atom 
are commonly treated to make two contributions to the hyperfine field Bhf 
on its nucleus [1]. The first is the contact Fermi field produced by s-electrons 
[2], which in the Gauss system of units is equal to 



^ = -y^|* s (0)| 2 <r. (1) 

The second is the dipole field created by the distribution of density of spin 
magnetic moment of electrons from the shells with non-zero orbital moment 
I [1,3] 

Bdip {tr ) = —jis J \*i(r)\ 2 3{V(T) r r ~ r2(T dr. (2) 
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In these formulas a nucleus is considered to be placed in the point of origin 
r = 0; fiB is the Bohr magneton; cr is the unit vector in the direction of the 
spin moment of electron; \1/ (r) - the electron wave function whose subscript 
in (1), (2) designates only the orbital moment of a shell. 

The two fields (1) and (2) are created by two parts of the same distribution 
of electronic spin moment density. Nevertheless, the contact Fermi field (1) 
apparently differs from the dipole field (2) in that it is dependent solely on 
the value of the wave function at r = whereas the latter is controlled by the 
distribution of spin density all over around the nucleus. And what is more, 
the Bp value is independent of cr, while the value of B dip depends on the a 
orientation relative to selected directions of the spin density distribution. That 
is why it seems appropriate to formulate the following questions herein. 

(i) Why does not the Fermi field depend on the distribution of spin density 
all over around the nucleus? 

(ii) Why is the multiple in Bp equal namely to 8n/37 

(iii) Would the Bp value change if the shape of \l/ s (r) changed whereas the 
\l/ s (0) value does not change? 

In what follows the answers for these questions are proposed based in part on 
the results related the field acting on a muon at interstice in solids [4]. 



1 The local contribution to magnetic induction 

In nonrelativistic quantum mechanics just as in classical electrodynamics, 
magnetic field of any magnetization distribution 1 

M(r) = -pL B \V(r)\ 2 <T (3) 

is defined by integral (2). If the point of observation r = lies inside the dis- 
tribution and the spin density does not vanish in this point, the under- integral 
expression diverges. Using standard techniques of calculating the above im- 
proper integral (see, for example, [3,5,4] ) and assuming that the magnetization 
changes over the space continuously, one can obtain expression 

^(0)=4 7 rM(0) + / r ^ ( r f {T)) dr, (4) 

1 In what follows we will not tell magnetization distribution from the distribution 
of spin magnetic moment density. 
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In the essence, this is simply another form of the standard formula B = 
AtiM+H that defines the magnetic induction inside the magnetization dis- 
tribution. The first term herein can be called a local because it, as well as 
the Fermi field, is proportional to the magnetization value at the point of ob- 
servation, though with another coefficient. The second term in (4) explicitly 
does not depend on M (0) and is determined only by the way of changing 
magnetization over ambient space 2 . 

Implicitly, however, the magnitude of integral contribution to (4) strongly 
depends on the M (0) value. If M (0) = 0, as is the case of spin density of 
electrons from atomic shells with / ^ (when this contribution coincides with 
(2)), M (r) first increases upon moving off a nucleus and then vanishes at r — > 
oo. As is easily seen, the regions of increasing and decreasing magnetization 
with r make contributions of different senses into integral in (4). As a result, 
the magnitude of hyperfine field is controlled mainly by the distance from a 
nucleus to the area where spin density is mostly localized, i.e., by the value 
of parameter (r~ 3 ). This designation in [1,6] was used for the radial part of 
integral (2), and with the designation for the angular integral part taken also 
after these authors, one can write the dipole contribution in the form 

(a ■ B dip (<r)) = -n B (r~ 3 ) (3 cos 2 9 - l) (5) 

Unlike is the case when M (0) 7^ and magnetization falls monotonically to 
zero with increasing r. Providing the distribution of magnetization is isotropic, 
M (r) — M (r) (s-electrons), it vanishes at r = where the first node of 
the radial s-function is located (for 1 s-electrons r^p — > 00). The integral in 
(4) over the region r > is equal to zero (see below), and in the range 
< r < r(p magnetization only lowers, so that this integral turns out rather 
large. Besides, because of the spherical symmetry, it does not depend on the 
form of function M (r) and always equals — -yAf(O). It is the sum of two 
terms in (4) that results in the Fermi field (1). 

What would change in the situation described if the part of spin density that 
does not vanish at r^0 is not of spherical symmetry? To answer this ques- 
tion it's enough to distinguish the two parts of the spin density distribution. 
The first of them is spherically symmetrical afresh and creates the same field 
87riVf(0)/3 at the point of observation. The second one has lower symmetry 
and vanishes in this point, and hence its contribution to the field is described 
by the formula (5). So the field at the point of observation now is not equal 
to Fermi field and, which is more, its value depends on M direction. 

2 In the classical electrodynamics, it is the latter contribution that, together with 
allowance for the discontinuity of magnetization on the body surface, determines 
the demagnetizing field. 
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For instance, in the case of the spin density distribution when the surfaces of 
constant magnetization have the shape of ellipsoids, a well-known result [5] is 
obtained from (4) 

B hf (r = 0, a) = 4tt [1 - N (a)] M (r = 0) , (6) 

where M (a) is the demagnetizing factor of an ellipsoid. In other words, the 
magnitude of integral in (4), still remaining proportional to M (0), turns out 
dependent also on the direction of spin moment, which is designated as pre- 
viously by symbol a. For the spherically symmetrical distribution of magne- 
tization this dependence disappears, M ia) = 1/3 so that An (1 — A/") = 8n/3 
and one comes again to the Fermi field (1). 

Let us discuss now how local expressions (1) and (6) arise from the integral in 
(4) . The matter is that a homogeneously magnetized layer that lies between 
two similar ellipsoidlike surfaces (in particular, spherical layer) does not pro- 
duce the field inside [5] . Representing the whole magnetization distribution in 
the form of an "onion" consisting of such closely adjoined layers, we can take 
off any external layers without changing field in the center 3 . As a result, the 
field in the center turns out coinciding with the field of an infinitely small ho- 
mogeneously magnetized ellipsoid, and the same conclusion follows from the 
common simple method of calculating Fermi field [3] as well. 

However, it does not mean that the Fermi field is produced solely by the 
infinitesimal central part of the spin density. It means only that for some 
spin density distributions of high symmetry (in particular spherical, cubic 
and ellipsoid-like symmetry) the sum of contributions to Bhf, produced by 
the outer parts of the distribution, becomes equal to zero. This is not true for 
the spin density distributions of lower symmetry, and to make matters worse, 
in this case one can not find unambiguously the part of the spin density that 
does not vanish in the point of observation, i.e. the part that should create 
the "contact" field. If, in addition to this, one takes into account the Bhf 
anisotropy, which is absent only in a spherically symmetrical case (1), to infer 
on the "locality" is hardly pertinent. 



2 Spin density anisotropy and Fermi field in crystals and molecules 

Could, however, one expect that in an atom there is a part of spin density that 
would lower monotonically and anisotropically upon moving off the nucleus? It 

3 It is just the reason for the part of spin density of s-electron from the region 
r > Tq 1 -* to not contribute to hyperfine field. 
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is impossible of course in an isolated atom, but it may be the case in a crystal 
of low-symmetry. Apparently, it is how, for example, the density of itinerant 
electrons in crystals with the lattice symmetry lower than cubic should behave. 
The simplest argument in favor of this quite evident statement can be obtained 
in the approximation of almost free electrons. Consider an electronic state 
= exp (ikr) in the rhombic r -type lattice. Let the perturbation potential 
has the form 



V(r) = 2V 



cos 



2nx 



+ cos 



+ COS 



2nz 



(7) 



If the V value is quite small comparing to the change of kinetic energy of 
an electron e(k) = 4p^, in the lowest order of the perturbation theory, six 



states 



kxi kyj k z ~\~ (X 



2n 



a 



±1 are 



kx ~\~ oi a , fcy , k z j , k x , ky -\- (x j , k z 
mixed to each state \k) with one and the same matrix element V. Hence, the 
charge p (r) and spin magnetic moment M (r) density of an electron acquire 
an inhomogeneous addendum proportional to 



2mV 



cos 



COS 



+ 



COS 



+ 



l) - (^) 2 (!) -(kyf (f) -(*.) 



(8) 



Hence, the distribution of spin density around the point of origin (the nucleus 
of one of the atoms constituent the lattice) is already anisotropic. In the vicin- 
ity of a nucleus, isolines of spin density turn out to have the form of ellipsoids 
elongated the stronger, the greater is the difference in lattice parameters a, b 
and c. 

Similar anisotropy should arise in the opposite limiting well, i.e. in the 

tight binding approximation, if to allow for the fact that the wave functions of 
neighbors overlap the stronger, the closer they are located. To make it evident, 
a study of the simplest case of a two-atom molecule is sufficient. Let us treat 
an ionized hydrogen molecule H2 in which atomic nuclei have the coordinates 
riand r-i, with the vector R = 7*2 — r\ being parallel to the z axis. The 
ground-state wave function is chosen in the form 



(r) 



1 

Q 



e -\r-n\/a _|_ e -\r-r 2 \/a 



(9) 



Q = yW 3 [1 + (1 + p + p. 2 /3) exp (-p)}. 



Here, p = R/a, a = a /Z is the radius of the Is orbit, which at R — > 00 
is equal to the Bohr radius a and decreases when the atoms come closer [7]. 
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This effect is taken into account via introduction of an effective nucleus charge 
Z whose value is conditioned at each R by the minimal average energy of the 
state (9). 

The spin magnetization distribution in this state is the sum of two contribu- 
tions: 



M{r) =Mi (r) + M 2 (r) , 

Ml (r) = U- 2 \ r -^\l a + e- 2 \ r - r '^ a ) , (10) 

M 2 (r) = ^g^ e" (|r " ri|+|r " ra|)/o . 
Q 

The magnetization distribution Af i (r) creates at every nucleus the magnetic 
field induction: 



B 



dip 



8tt 1 + e~ 2p 
3 (nR) R - \xR 2 



Bi = — -Hb<t — 7^ — + B dip , (11) 
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This result can easily be obtained after simple reasoning similar to the above 
with account for the fact that a homogeneously magnetized sphere creates 
around the same field as a point dipole placed in the center of the sphere and 
having magnetic moment equal to that of the sphere. 

The formula (11) looks like the standard expression for calculation of the field 
at an atomic nucleus in a crystal, which includes the Fermi contribution and 
the dipole fields created by the neighboring magnetic atoms. However, in this 
case, the homogeneous and equal to aexp (— 2p) /Q 2 part of the spin density 
distribution, which is concentrated inside the sphere with a radius R and the 
center located at the neighboring site, contributes solely to the Fermi field. The 
dipole field is given rise to only by a heterogeneous part of this distribution 
so that the magnetic moment entering into B^ip (11) is equal 



A-kiiec f R 



/ x 2 (e- 2x/a - e- 2p ) dx (12) 
Cj Jo ' 

1 1 - (1 + 2p + 2p 2 + 4p 3 /3) exp (-2p) 

~2 <T/ii? i + (i + p + p2 / / 3)exp( _ p) 

As a result, with decreasing R, the anisotropic contribution B dip reaches its 
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maximum and vanishes at R — » 0, rather then increases proportionally to 
1/R 3 . 

The second contribution to the magnetization M 2 (r) arises because of the 
interference of the atomic states in (9). The surfaces M 2 (r) = const are 
the ellipsoids of revolution, with the eccentricity tending to zero at |t* — t*i| H- 
\r — 7*2 1 — > R. Again, with the considerations similar to the above-employed, 
it is easy to show that the contribution of M 2 into the magnetic field induction 
has the same value and direction at the nuclei and at every point of the segment 
that connects them. Just as the Fermi field, this contribution is proportional 
to the magnitude M 2 (r) , which is constant at every point of the segment and 
equal to —2p B aexp (— p) /Q 2 . Therefore, the field induction related to M 2 is 
determined by the formula similar to (6): 



M(a\\z) = -l{p)-2^ (13) 



B 2 (a) = -Sitp B a [1 - M (a)} exp (-p) /Q 2 , 

roo / x _|_ 9 s 

l(p)=p J e~ px x (x + 1) (x + 2) In ( —— ) dx. 



The dependences of B\ and B 2 on R/ao at er\\z and a^Lz are shown in the 
Figure 1 and corresponding dependences of the anisotropic contributions to B\ 

and B 2 , i.e., B dip (cr) and B Fan (cr) =B 2 (a) - B 2 (ff\\z ^ +2B 2 (a±.zj /3, 

are displayed in the first insert. The total anisotropy value makes up about 
10% \B 1 + B 2 \, and the Fermi-field anisotropy turns out only slightly lower 
than the dipole field anisotropy, as is seen from the second insert in the same 
figure. 

Thus, in the H 2 molecule, with decreasing spacing, not only the hyperfine field 
at the nuclei as such increases (chemical shift) but its anisotropy as well. In 
this case, about half the anisotropy contribution is ascribed to the anisotropic 
part of the Fermi field related to hybridization of the atomic orbitals. The 
same effect should evidently arise also in crystals in which the distances from 
a magnetic ion to its closest neighbors in different directions strongly vary, 
which leads to various overlapping of their outer s-shells. 

In general, the value of local contribution to the induction (4) is controlled by 
the form of dependence div (M) on coordinates. If the spin part of the wave 
function does not depend on r, i.e. M (r) is directed parallel to ^-axe in any 
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1 2 3 

R/a Q 



Fig. 1. The dependences on R/aQ of two contributions into the hyperfine field at 
nuclei of the molecule - B\ (11) and B2 (13) for different direction of the spin 
a. In the first insert, the dependence of the anisotropic parts of these contributions 

- Bdip(a) and Bp an (a) =i?2(c) — 167r//exp (— p) / (3Q 2 )on R/clq are shown, while in 
the second insert, the ratio Bpan/Bdip independent of a. 

point, div (M) has the form 

.. ,„., „ dM n ldM . n 

div(M(r)) = ——cos 9 — smS. (14) 

or r o9 

The first term in this formula having been integrated on r in (4) always gives 

— ^-M (0). The result of the integration of the second term in common case 
is the functional of M (r). In particular, if M (r) = M. (r) $ (9, </?), as is the 
case for atomic states with / 7^ ; the value of this functional depends on 
(r~ 3 ) just as in (5). For the result of this integration to be, alike the Fermi 
field (6), proportional solely to M (0), the M (r) dependence should have the 
special form 

M(r)=M(r-$(e,(p)) (15) 
where M. and $ are any sufficiently smooth functions. Taking into account 
that in this case 

ldM _dM 1 9$ 
we obtain replacing (14) and (16) into (4): 



B(0) = M(0) 



sin# <9$ 

cos 9 — 



$ (9, cp) 89 



sm9d9dtp\. (17) 
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The atomic s-state is a trivial (the angular dependence is absent) particular 
case of spin density distribution (15). The other examples of distributions of 
this type are ellipsoid-like distribution or a distribution with cubic symmetry. 

In order that the spin density distribution for an isolated atom turns out a 
function of the type (15), its ground state should be mixed to by the states 
of the continuous spectrum. Such mixing up actually arises when an atom 
is located in an electrical field, in particular, the crystalline electrical field. 
However, calculation of the corresponding contribution to the Fermi field is a 
serious mathematical problem, which requires special treatment. 



3 Conclusions 

Returning to the questions posed at the beginning of the paper, we now can 
formulate the following conclusions. 

(i) The contact contribution to the induction at internal point of any magneti- 
zation distribution virtually results from the summing up of the fields created 
by all elements of this distribution. 

(ii) The value of the proportionality coefficient between the induction and 
magnetic moment density at a point of observation is controlled by the sym- 
metry of the magnetization distribution. This coefficient is isotropic and equal 
to 8n/3 only for distributions with a spherical or cubical symmetry. 

(iii) In molecules and crystals of a lower symmetry owing to the hybridization 
of the states of neighboring atoms, the symmetry of the s-electron density 
distribution around the nucleus decreases. This results not only in the propor- 
tionality coefficient between the induction and magnetic-moment density at 
the point of observation changes, but starts depend on the moment direction. 

It should be underlined also, that the main result of this paper reduces to the 
substantiation of a very simple statement: individualization of the Fermi field 
from the total hyperfine field becomes empty of meaning for low-symmetry 
crystals and molecules. The calculation of the field created by a spin density 
distribution of low symmetry should be performed by the formula (4) or some 
equivalent to it. The standard presentation of hyperfine field as the sum of the 
Fermi (1), dipole and Lorentz fields and demagnetizing field of the specimen 
results in masking the role of the spin density inhomogeniouties, which may 
be the source of some faults. 

Allowance for this circumstance may be quite important when analyzing the 
results of the hyperfine field measurements. The case in hand may be both the 
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fields acting on a probe particle (neutron, muon and others) in the interstices 
[4], and the field on an atomic nucleus in a crystal or molecule. In all these 
cases, with a proper account for the contribution of collective electrons and/or 
distortions of the atomic s-shells, additional anisotropic contributions to the 
Fermi field can emerge. In particular, one could not exclude that the high (up 
to 3 -J- 4 T) anisotropy of the hyperfine field at the nuclei of iron atoms that 
constitute the dumb-bells in the compounds R^Feyj [8,9,10] is related to this 
effect rather then to unquenching of orbital moments. 
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